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ABSTRACT 

Ten basic stops are established fur an analog mcth- 
oii that measures runt ml sxstcm response parameter*. 
An example .shows how these steps were uBed on a speed 
control portion of an uuxllinr) |s>wer unit The equations 
and i alculntions necessary to describe this subsystem 
are given. The mechanization schematic and simulation 
diagram for obtaining the measured response parameters 
of the control system using an analog cireiut are ex- 
plained Methods for investigating the various effects of 
the control parameters are described. It is concluded 
that the optimum system should Ik- underdamped enough 
to be slightly oscillatory luring transients 

INTRODUCTION 

Control system response ps. amidol's can be ob- 
tained by various techniques, some of which a ro listed 
below: 

1 Check out of subsystems with analytical system 
studies. 

I! Tes* equipment and phvsicnl duplication of the 
controlled process 

3. Digital simulation of the controlled process. 

I. Analog simulation of the controlled process. 

Method one Is simple and does not require complex 
apparatus nr setup time. Its mujor disadvantage is 
customer dissatisfaction System startup is quite d'i- 
ficult when no measurements have been made on the en- 
tire system prior to setup in the operating site Much 
time is required to "debug" the controlled process and 
get it to meet the customer specifications 

Method two is complex and requires working with 
complex apparatus usually not available. It is prohib- 
itively expensive to set up a complete system test for 
all possible combinations of processes and equipment. 

For these reasons many manufactures make measure- 
ments on ns much of the controlled process as possible 
and hope that no insurmountable problem occur at the 


job site 

Method three Is simple and suitable for solution of 
nontrivial mnthmatie-'l equations Several software 
packages are currently available with which controlled 
processes can be analyzed iRef. 1) Its major di sad van 
tugc is that a digital computer is required to (icrform 
the nnulysis which may very well not be available to 
many designers. If you have aceess to a digital com- 
puter, bv all means lo !< further into this method of 
analyzing controlled processes. 

Method four, the subject of this tutorial type pa|ici , 
is generally a good compromise betueen the other meth- 
ods. Analog simulation is not complicated after one 
becomes familial with how it Is done. Circuit compo- 
nents are available either at a reasonable cost or by 
renting. 

There are two basic approaches to electronic analog 
simulations. The conventional approach is to use a 
general purpose analog cumputer. In this ease the 
equations would be written in a form that contains only 
Integrators, summers and gain elements. Direct iden- 
tity of physical parameters and transfer functions are 
often lost in the changes necessary to construct a sim- 
ulation that contains only these elements. The second 
abroach presented here shows how to produce a sim- 
ulation with an electronic circuit constructed with op- 
erational amplifiers for minimum cost. This approach 
could be us-d even If your faellity is very limited in 
available funds. 

There are ten basic steps in enveloping analog slm- 
ilation for measuring control system response. The 
steps are as follows: 

1. study the control process 

2. Draw a sy item diugram of the process 

3. Obtain transfer functions for each block 

t. Kstimatc the process constants 

5. Calculate netw< k simulation constants 

G. Arrange equations to avoid differentiators 
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7. Derive the network simulation schematic 

8. Dele t mine magnitude and time scaling foe the 

variables 

a. build the s|ie< l.ill/cd simulation circuit 

10. obtain simulated system response data 

DESIGN EXAMPLE 

As on example o( how the steps produce solutions 
to practical problems consider the speed control portion 
of an auxiliary power uni 1 API/). Al’l' control systems 
generally sense output voltage, current, phase and fre- 
quency. These variables art then com|Mired with fixed 
references to determine the power systems electrical 
performance. In our sample portion of these larger 
systems, a throttle valve limply controls the rotating 
shaft speed of the alternator by regulating the rate of 
fuel consumption. Fig. 1 shows the speed control sub- 
system dingrum. From a study of the control processes 
Involved in this portion of the hardware or from theii 
s i ice if hat ions, one con determine how the various blocks 
are connected together and prepare a diagram of the 
process. This simplified speed control shows u turbine- 
driven alternator whose output ac waveform is controlled 
by a combination of three feedback loops - phase, fre- 
quency and throttle valve position. 

The transfer function diagram is derived from the 
subsystem diagram and from knowledge of the compo- 
nent dynamics. For example, the self- saturating mag- 
netic amplifier in the subsystem diagram, block number 
3, is known to have an approximate iransfer function of 
the form shown In Eq. 1 (see refs. 2, 3, und -1). 

E 0 (s)/I c ts> » Kj/dj s + l) V/mA (1) 

The constants K. ( and r,, are obtained from specifica- 
tions, calculated or measured experimentally. Consid- 
erable work has been done to develop appropriate trans- 
fer function for various devices. Itefs. 5, 6, 7, and 8 
should give the reader a starting point from which to 
obtain transfer functions for most control devices. One 
can, however, always go back to the original differential 
equations and derive whatever is required. 

Fig. 2 shows the transfer function diagram for the 
speed control process. An index is used to identify the 
control process components by numbei . For example, 

Kj represents the gain of the 30 Hz phase network and 
K., represents tha gain of the pulse time modulator. 

The fifth step in the technique is to obtain the proc- 
ess constants. Again, these constants may be obtained 
from the control process by measurements on the com- 


ponents, from calculations or from the specifications. 

For our example, consider transfer function nine which 
represents the subsystem circuit shown in Fig. 3. This 
lead- lag network has the transfer function 

E 2 is)/Fj(b) • Kjjil ♦rjjS) '(l ♦TyS) (2) 

Source and loud Impedance have been neglected on the 
assumption that It. « Itj « R^. The equations and 
calculations necessary to define Ky, r y, and r. are 
also derived in this figure. The network constants are 
given by the initial specifications to be Hj ■ 300 K;’, 
Cj“0.18uF and It,, ■ too K'l. Using the equations 
derived in Fig. 3, calculations for the lead-lag network 
constant* result In Ky “ 2.5*10"* V/V, 5. 4xlo"*s, 
and Ty * 1.3x10*" *. Later on In this technique we may 
want to vary these constants, if it turns out that the po- 
sition loop is a cause for Improper performance. The 
remaining process constants are shown In Table 1 and 
were obtained in a similar manner. 

BUILUfNw THE SIMULATION NETWORK 

The analog simulation schematic diagram can be 
prepared from the transfer function diagram by com- 
bining active networks. Fig. 4 shows the In, ut-nutput 
voltage relations for six basic circuits that could be ap- 
plied to some problem. Additional transfer functions 
simulation networks can be found in Refs. P to 12. 

Let us explain how transfer functions would be 
mechanized by first modeling the servomotor and gear 
train. These components have a transfer function as 
given in Kq. 3. 

Oj/Ij « K |/s(l + T,, s) (3) 

Combining transfer functions nun her 3 and 4 from 
Fig. 4 gives the required simul'.don as shown in Fig. 5. 
The resistor and capacitor values are obtained using the 
equations given in Fig. 4 and the constants from Table II. 
In this manner all simulation constants are calculated 
for the complete schematic. 

In preparing this diagram one should also give con- 
sideration to the variables of the control process. The 
magnitude and frequency range of these variables are 
scaled to fit within the *15 V range and the 0 to 100 Hz 
frequency range allowed bv the |x-n recorders. The Al’U 
performance is slow enough that it was studied in real 
time and the variables were all limited and scaled as 
required to correspond with physical control boundary 
conditions. Each monitored variable will require an 
additional channel of recording equipment so some care 
must be used in selecting the variables which require 



monitoring. For our example, live variable* defined 
Ix-low required monitoi mu. 

Monitored Process Variable* 
b> rnbol Definition 

u j Valve position, degrees 

I*. Decomposition ehamber pressure, psig 

T, Loud toiquc, lb- ft 

i 8 Phase shift, radians 

ifeg/dt Frequency, ruiUons/ second 

The frequency network (transfer function thirteen) 
suggests that . H be differentiated, UlPerentlators ure 
normally avoided in simulation problems as they are 
peaking circuit* tending to overdrive (suturatc) ampli- 
fiers. Fortunately, all stable linear systems can be 
simulated without differentiators, see lief. 13. Non- 
linear simulations which demand differentiation can be 
approximated by (he circuit in Fig. I if phase load and 
gain limitations at high frequencies are aciept le 

Since our simple examp.e is linear, we avoid this 
complication by rearranging the circuit. We are able 
to obtain d s /dt directly by combining K-, and K 8 with 
a simple lag ns shown in Fig. 2. Integrating this signal 
gives us ,/g. 

Analog simulation has several known sources of 
error (Ref. 15). Methods described in the literature can 
be used to Identify how accurate the simulation will be - 
typically around 90 percent. Even with accuracy limita- 
tions analog simulation is an excellent method for pre- 
dicting where problems might occur ami for planning 
methods by which corrective action may be implemented. 

Its low cost is also a key advantage. We could buy 
the simulation network coni|ionents for n reasonable cost 
estimated below. Notice that these parts could easily 
be purchased for less than $250. 00. To make permanent 
records for each condition studied in thit problem, we 
could rent a five channel recorder for a few months for 
less than $500.00. This whole problem could be run for 
less than $ 1000.00 - a small price to pay for all that is 
being learned. 




Number 

Cost 

1. 

bummer Amplifiers 

5 

$50 

2. 

Integrator Amplifiers 

8 

80 

3. 

Potentiometers 

5 

25 

V, 

Power Supplies (dual) 

1 

80 

235 


ANALYZING THE DATA 

By applying a step-change in load toiquc to the sim- 
ulation network through an amplifier we cun observe how- 
control variables will respond for various conditions In 
the network simulation. Fig. ti shows transients of sev- 
eral variables when n 40 hp step load is applied to the 
s|iecd control. The top graph in this figure plots valve 
IKisition as a function of time. The valve would hit the 
position limit at 40° hi about 0. i s when the load is 
switched on. Thus protection against damage due to 
rapid slamming will be required. After a 1 second 
oscillatory condition, the vnlvc is open about 32° w ithin 
the design operating range. Chamber pressure, I’., 
reaches a 320 psig peak during transient periuds. Re- 
liable design practice requires that the decomposition 
chamber be designed with a safety margin of 5 for these 
|ieak stress conditions. The steady state chamber pres- 
sure would lx- about 230 |>slg. Graph I shows that the 
|x-ak frequency error, d, 8 /dt, is about 3.75 rad/ sec. 
and occurs at about 0.25 sec after the load Is switched 
on or off. The frequency error goes to zero lr. about 
2 sec. The maximum phase error, ,, 8 , is about 17 rad- 
ians as can be seen in graph 5. The speed control with 
the assigned constants does not exceed the allowable 
phase error, approximately 25 radians, aid therefore, 
would not drop out of synchronization for t ils loud 
change. 

Individual run can be made to determine the effects 
of varying gains, process constants, environmental 
factors or aging on the speed control performance. 
Variations of any of the netw-ork constants in the simu- 
lation can be accomplished by inserting new values for 
tile resistors and capacitors us required. I’otentlom- 
eters make tnese variations much easier If your circuit 
design allows for their use. For each set of conditions 
a recording like Fig. 6 is mniie and evaluated. The 
response shown In Fig. 6 Is underdamped enough to be 
slightly oscillatory. Further Investigation showed, 
however, that this was about optimum for this speed 
control with its inherent decomposition lags, see 
Ref. 16. 

CONCLUDING REMARKS 

The analog method of system design was accom- 
plished through ten basic steps. This design process 
was illustrated by analyzing a speed control portion of 
an auxiliary power unit. The equations and calculations 
nec essary to desc rlbe the subsystem are given. The 
simulation schematic and diagram for obtaining meas- 
urements of the control system response using an analog 
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circuit* urc explained Once tic simulated tontrul 
process is setup in network fnrm, Individual runk can 
easily be made to evaluate I be eflectk of varying gains, 
protest tonstanls, envlrcnmcntal fat tors or aging fav>- 
tors on the control •yktent performance. Variation* ol 
any of the network constants In the kimulalion can be 
a« omplUhcd by inserting m u values for the resistors 
ami capacitors. For each set of conditions a lime re- 
sponse rciording Is made and evaluated. The system 
should be underdumped enough lo be kltghtly osiillatorv 
during transients to optimize the control process. 

REFERENCE'S 

1. F. II. hpeci.'tart and W. L. Green, A Guide to I sing 
CS.MP-T 1 : ronttnuous bystem Modeling I'rogram- 

A [irog.am for simulating physical systems . 
Englewood Cliffs, N.J.: I’ren tier -Hall. 1976. 

2. A. krigman, "Circuits for Transfer Functions," 
lastrum, Conlml Syst . . Vol. 42. pp. 67-09, 1969. 

5. J. Ci. Truxal, Automatic Feedbm k Conttol System 
Synthcsb , New York: McGraw-Hill, 1955. 

I. W. It. Ahrcndl, Servomechanism I ractices . New 
York McGraw-Hill, 1954. 

5. II. K. Baffin and It. Muffin, "Approximating Trans- 
fer Functions From Frequency ltrs|»>nac Data," 
Instrum. Control Syst. . vol. 58, pp. 157-144, 1969. 

6. D. J. Mmmons, "Finding Transfer Functions From 
Pulse Response," Electron. Design , vol. 12, 

pp. 72-77, 1964. 

7. W. F. Cowherd II and T W. Cadman, "Numerical 
Prediction of Transfer Functions," I nstrum. Con - 
trol .Mad . . vol. 42. pp. 109-1 15, 1969. 

8. A. .Suzuki, K. Cotoh and T. Tanaka, "Obtaining 
Frequency Characteristics from Transient 
Resimnse." Instrum. Control Mat ., vol. 10, 
pp. 114-117, 1967. 

9. F'. F. Nixon. Principles »f Automatic Controls . 
Englewood Cliffs, N. J.: Prentice- II all Electrical 
Engineering Series, 1955. 

10. W. J. Karplus, Analog Simulation; Solution of 
Fluid Problems , New York: McGraw-Hill, 1958. 

11. Applications Manual for Computing Amplifiers for 
Modeling, Measuring, Manl|iulating and Much Fisc , 
Boston: Philbrlck Researchers, Inc., 1966. 

12. Handliooh and Cataloged Operational Amplifiers . 
Tucson. Ariz.: Burr- Brown Research Corp. , 

1969. 


15. J J. D'Azxa an<l C. Iloupls, l.lneai Control s\ stem 
Analysis anil Design Conventional and Modern . 

New York: McGraw-Hill, 1975 

II. V. It. Lalll, "Servomechanism Design Techniques 
and Abdications Aerospace Problems," NASA I'M 
X-6H25I, 1975 

15. R Stain, "Operational Amplifiers, Parti- 

Principles of Operation and Analysis of Errors," 
Electromechanical Design , vol. 9, Sep. 1965. 

16 C. W. Fetheroff, "Analysis :uid Computer btud\ ol 
tlu' (1-26 At I Speed Control System," TRW I H- 
2910, 1955. 



TAHl.K I - 1’ROCESS CONSTANTS 


Symbol 

Nominal magnitude 

Mange 

I’nlta 

K 1 

k 2 

1 

0 to 2 

Radian/ radian 

Reference, 1.43 
I’hane. 1.43 

0 to 1 . 91 

mA/ radian 

K 3 

IH>sition. 20. 9 
Frequency, 7.1 
I’haue, 0.028 


V/mA 

»3 

0.01 


Sac 

K 4 

Motor, 0. 0015 
Gear train. 0.021 


Dyne- cm/ rm a V 
Kadi an/ radian 

»4 

200 


Dyne- cm/ radian/ aec 

J 4 

3.6 


Gm- cm 2 

k 5 

7.5 


1W degree 

t 5 

K 6 

75 

50 to 100 

Msec 

0.0114 

0 01 to 0.0128 

lb- ft/ pal 

i/d 7 

3580 

1920 to 5240 

RatV aec/ lb- ft 

t l 

1.375 

0 to 2.75 

lb- ft 
Mug- ft 3 

J 7 

0.00235 

0.001 to 0.0037 

k 8 

0.314 


Radian/ radian 

K 9 

0.25 


V/V 

t 9 

0.054 


Sec 

t 9 

0.0135 


See 

K 10 

K n 

1.81 

1.08 to 1.81 

mV rms V 

0.20 


Rms V/dcgrec 

Kj£ 

0. 154 


V/V 

t 12 

0.25 


Sec 

t 12 

0.01 


Sec 

K 13 

0.995 

0 to 1 99 

mA-scc/ radian 

0. 125 



Radian/ radian 
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TABU: II - COMM TKIl CONKTAN'K 



Nominal magnltuilv 

Pol Range 

Holla 

b hematic parta 

K t 

K 3 

1 

0 to 0.23 

Haitian/ raillan 

P. 

Reference, 1.49 
Iliaee. 1.43 

0 to 0.76 

m A, radian 
V/mA 

1*2 H4/K2 R4/IU 

•4 

Modulator, 3.3 
Poaltton, 3.3 
Frequency, 37.2 


V/mA 

K3/H3K3C2 R3/MR9C3 
RA/R7R9C2 

Tj 

0.01 


tec 

Cl Ra 

*4 

9.23 


D»|rw/ acc 
Kma V 

R11/R10 


o.ois 


Mr 

R11C3 

K i 

T.» 


Pal/ degree 

Rl»/R17 


0.075 


tec 

R18C6 

K 6 

0.011 

0.10 to 0. 16 

lb- ft/ pal 

Ps 

k 7 k„ 

3580 

0.30 to 0 63 

Kad/irc 
lb- f 1 

H21/R11. .I21/R20 

» T 

13.23 

0.30 to 0.83 

tec 

R21C7 


0 453 


mA rmi V 

Rid/ HI 5 

r 9 

0 054 


tec 

RISC! 


0.013 


tei' 

R14C3 

K„ 

0.20 


Kma V/degree 

R14/R1S 

Kj2 

0.038 


mA/roA 

H2T/R25 

r \2 

0.23 


tec 

R25C0 

r 12 

K 13 

K lt 

0.01 


tec 

R27C10 

005 

0 to 0 00 

mA/mA 

R24/R23 

0.125 


Radian/ radian 

R31/R30 
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APU speed control subsystem diagram. 
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Figure 3. - Lead-lag network circuit constants. 
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CHOOSE RIO - 100 KO 
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Rll 

Rll * 9.25 x 100 K - 925 KO 
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U* 
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C3 1.8x10*2 sec 
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t 4 ■ 1. 8xl0' 2 sec 
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•a 0194x10 6 F 


Figure 5. * Network and constants for transfer function 4. 
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